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Abstract 

Nuclear muon capture on the proton and ^He is considered both within and beyond the 
Standard Model in terms as general as is possible. Explicit and precise analytic expressions 
for all possible observables are given, assuming only a Dirac neutrino in the limit of vanishing 
mass. These results allow both for precision tests of the Standard Model and new physics, as 
well as for the assessment of the potential physics reach of experiments designed to measure 
specific observables. Using these expressions, stringent constraints can already be inferred 
from a recent precision measurement of the statistical capture rate on ^He. Likewise, similar 
constraints should follow the completion of a precision measurement in progress of the singlet 
capture rate on the proton. 
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1 Introduction 



Over the last fifteen years, muon physics has regained its rightful place in particle physics, rang- 
ing from intermediate energies well into the high energy frontier in the foreseeable future with 
the advent of muon colliders. Given the availability of intense muon beams at different labo- 
ratories, as well as new and much efficient experimental and detector techniques, intermediate 
energy muon physics has moved into the realm of precision studies of the Standard Model, with 
the hope of possibly unravelling some tell-tale sign for the physics which must lie beyond it 
in ways complementary to present day high energy experiments at the colliders. Such studies 
include both purely leptonic as well as semi-leptonic electroweak processes, the much studied 
field of nuclear muon capture Q belonging to the latter class. 

When the aim of such electroweak processes in nuclei is to study particle physics issues, the 
uncertainties inherent to nuclear structure modeling have to be disposed of to the largest extent 
possible, leaving essentially only the lightest of nuclei available, beginning with the proton and 
next the stable 3-nucleon bound state, ^He, since indeed the quantum wave functions for the 
latter system are by now very well understood |2|, ^. An additional bonus in the case of these 
two systems is that they both provide the unique occurrence of a (mother-daughter) pair of 
nuclei which define a spin 1/2 isospin doublet, namely the proton and neutron, and ^He and ^H, 
respectively. Hence in the limit of exact isospin symmetry, some of the various phenomenological 
nuclear form factors which parametrize matrix elements of these states are related to one another, 
while for spin 1/2 states, the numbers of these form factors remains small. This is how through 
CVC, the experimental determination of electromagnetic form factors from electron scattering 
may be translated into knowledge of the related nuclear form factors for the charged electroweak 
vector current. In other words, hiding our ignorance of the microscopic dynamics at the quark 
level into a phenomenological parametrization in terms of form factors, it remains possible to 
consider predictions of observables which do not require models for nuclear structure. In the 
field of nuclear muon capture, this description corresponds to the so-called "elementary particle 
model" approach to nuclear electroweak processes^, which will be used in this paper. 

One of the main motivations for studying nuclear muon capture, especially in light nuclei^, 
has always been to measure the induced pseudoscalar nucleon form factor gp{q'^), certainly the 
least well-known of all non vanishing nucleon form factors with a combined uncertainty which 
has stood at 22% for the last twenty years[^. The urgency of this specific issue has recently 
become more pressing, mainly for two reasons. On the one hand, based on the fundamental chiral 
symmetries which, even though dynamically broken, survive non linearly in the nonperturbative 
low energy regime of the theory for the strong interactions among quarks, namely quantum 
chromodynamics, definite predictions for gp with a precision of a few percents or better have 
been achieved||6|, 0, The experimental confirmation of the expected value is thus a crucial 

low energy test of our basic concepts for the theory of the strong interactions. On the other 
hand, a new 8% precise measurement of gp has been completed in the intervening years, using 
the rare process of radiative muon capture on hydrogenjlO, The unsettling fact is that 
the value obtained differs significantly (by a 4.2 a deviation) from the expected value. No 
theoretical explanation for this discrepancy having been found so far (see for example Ref.|^ 
for references), the issue thus remains open and this difficult experiment shrouded with some 
feeling of uncertainty. Indeed, the general understanding of the dynamical breaking of the chiral 
symmetries of QCD has so far never been found lacking in essentially any other low energy 
process. 



^In heavier nuclei, the same axial electroweak probe enables to address the issue of hadronic coupling and mass 
renormalization in the nuclear medium. 
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During the same period, another precision measurement of nuclear muon capture on ^He 
has also been completed 1 12]. In spite of its remarkable precision of 0.4% for the statistical 
capture rate at Ag**^ = 1496 it 4 s~^, whose value is found to be in perfect agreement with the 
theoretical prediction Xl\^^^^ = 1496 ± 21 s-^g, |, the ensuing value for gp, though at exactly 
the expected level, is still precise to only 19%]l3[. One should add however that in terms of the 
corresponding pseudoscalar nuclear form factor Fp{q^) for '^He, the same experimental result 
translates [H, |l|, |l5|, |l6|, [l^l into a 13% precise test of the prediction based on PCAC, the 
symmetry which together with CVC has historically been the precusor for the chiral symmetries 
of QCD. Even though these results are still a long way off the precision reached by the latest 
theoretical analyses j^, @, H, they stand as a clear confirmation of the basic concepts involved 
in the chiral symmetry aspects of the problem. 

Independently of this specific situation with respect to the value for gp, if one is willing to 
use the theoretically expected number, the latest ordinary muon capture experiment on '^He has 
reached such a level of precision that other tests of the Standard Model (SM) for the electroweak 
interactions become possible [||, l^, 0, some of which prove to be quite stringent for possible 
new physics beyond the SM. Moreover, given the above issues surrounding the value of gp, as 
well as the precision of its theoretical prediction, a new effort has been launched fl^, |l^ ] in order 
to measure the singlet rate of ordinary muon capture on hydrogen in a gas target (to avoid 
the complications due to molecular binding effects), hopefully to a precision of 0.5% to 1%Q. 
Here again, beyond the initial aim towards the value of gp, this level of precision should also 
enable tests of the SM. However, in order to infer a value for gp from any experiment, analytic 
expressions for observables whose numerical evaluation is to the required standard of precision 
should be available. This is the purpose of the present work, much in continuation of that of 
Refs.|,i. 

This paper considers the ordinary nuclear muon capture process on a spin 1/2 isospin 
doublet in terms as general as is possible. Not only are all possible effects existing within the 
SM included in the present analysis, but any possible contributions which may appear beyond 
the SM due to new interactions are accounted for as well through an effective four-fermi quark- 
lepton interaction which describes all possible scalar, pseudoscalar, vector, axial and tensor 
couplings. Indeed, given the momentum transfer involved, much less than any of the mass scales 
characterizing such possible interactions, an effective four-fermi interaction at the quark-lepton 
level is perfectly justified. 

Previous analyses over the years, beginning with Ref. [p^] pointing out the existence of 
the hyperfine effect in the capture rates because of maximal {V — A) parity violation^ have 
all only included vector and axial interactions through the usual {V — A) charged electroweak 
coupling, and only some of these works have considered the possibility of so-called second- 
class currents. A handful of observables have been computed, and then not always in analytic 
form[Q, |2^, 21, |2^, |2^, |2|, In contradistinction, the present work does not involve 



any non relativistic expansion and allows for all possible contributions, including CP- or T-odd 
effects through complex form factors and couplings coefficients in the effective interaction. The 
only approximations made are, on the one hand, that the neutrino mass is taken to vanish and 
that the muon leptonic flavour is conserved, and on the other hand, that nuclear recoil effects 
associated to the scalar, pseudoscalar and tensor interactions beyond the SM are also ignored, 
since the latter couplings are necessarily small on their own and would in turn multiply small 



^The principle of the experiment lies in the comparison between the muon disappearance rates, through the 
usual electron decay mode, for both positive and negative muons, only the latter being subject to nuclear capture. 
Since this experiment does not measure the neutron recoil energy distribution, it has no handle on the neutrino 
energy spectrum nor on its mass. 

^The same effect arises of course also for a pure (V^ + A) coupling. 
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nuclear recoil contributions. Moreover, all possible observables become available through our 
results, hopefully making it easier to assess their sensitivity to whatever parameter or physical 
input a specific experiment using ordinary nuclear muon capture on these nuclei would wish to 
address. This point will explicitly be illustrated with some of the obtained observables, in the 
context of the two experiments mentioned previously. 

Most of these observables require correlation measurements using polarized states for the 
initial muonic atom, while another subset requires to measure also the polarization states of 
the outgoing nucleus and/or neutrino. These are certainly experimental challenges bordering 
on the impossible for certain of these polarization observables, but some experimentalists take 
up the task. For example, there exists a first-generation experiment! 25 1 which measures the 
vector analyzing power of the outgoing triton in polarized muon capture on ^He. Even though 
the preliminary results are yet not precise enough to be of use in a theoretical analysis, they 
certainly demonstrate that this specific difficult challenge can be met. 

More specifically, the physics reason for these difficulties is that during its atomic cascade 
down to the muonic atom ground state, the muon, even if initially polarized, suffers depolariza- 
tion effects to a great extent, leaving over only a small fraction of its initial polarization|l|, 
p7| , p^ j. Once in the ground state, the degree of polarization may be increased again by external 
means! 



29, 25], but not to any large degree and less than theoretically anticipated [|30| . In fact. 



both the initial nucleus and muon need to be polarized in order to end up with a muonic atom 
ground state polarized to any degree] 31]. Hence, because of these atomic physics issues, even 
though the expressions for all observables are now available in analytic form, any experiment 
which consists not only in a rate measurement is in essence extremely difficult to perform when 
aiming towards the demands of great precision^. 

The outline of the paper is as follows. In the next section, the general parametrization of 
the capture amplitude is described, and the ensuing expressions for observables explicitly given. 
In Section 3, the situation for ^He is then considered in detail, first within the SM, and then 
beyond the SM, each time by using the experimental result of Ref.]12] for the statistical capture 
rate and also by illustrating the potential physics reach of some other final state distribution 
which does not entail a correlation observable except for an initially polarized muonic atom. 
In Section 4, the same considerations are applied again to the case of capture on the proton. 



based on the optimistic aim of a 0.5% precision in the result for the singlet capture rate]]!^, 14]. 
Section 6 then translates some of the limits for physics beyond the SM obtained from the 
previous considerations, into limits for parameters of some specific models for such physics. The 
Conclusions end our discussion, while further information relevant to the analysis is provided in 
three separate Appendices. 



2 Muon Capture Observables 
2.1 Kinematics 

Our notations for kinematics are as follows. Let m^, Mi and M2 be the masses of the muon, 
and of the initial and final nuclei, respectively, with M = (Mi -|- M2)/2 the mean value of the 

''in this respect, the recent proposal of Ref. appears to be totally unrealistic, the more so since some of the 
possible contributions which have not been included in that analysis, such as recoil order effects, could also lead 
to contributions to the massive neutrino polarization states. Such effects are all accounted for in our analysis in 
the limit of a massless neutrino, including the possibility of T-violating couplings and currents other than V and 
A. Note also that due to helicity constraints, transverse contributions to a massive neutrino polarization state 
produced in muon capture through only V and A couplings are necessarily suppressed by the ratio of the neutrino 
mass to its energy, and are thus unobservable. 



3 



latter two. As mentioned previously, the capture process is considered in the limit of a vanishing 
neutrino mass and no leptonic flavour mixing. Initially, the muon and nucleus form a muonic 
atomic bound state at rest, whose total rest-mass is denoted y/s, which differs from + Mi 
by the binding energy — (aZ)^/i/2, fi being the reduced mass of the bound state (expressions 
are given in units such that c = 1 and h = 1 throughout). Let then phe the momentum of the 
outgoing nucleus {—p is thus that of the outgoing neutrino), io the energy of that nucleus and v 
the energy of the neutrino, such that we have 

u=\p\ , V^=u + u , -=^^ , '" = ^7^- 

With respect to polarization states, let si denote the normalized polarization vector of the 
initial spin 1/2 nucleus, as measured in its rest-frame, and §2 that of the final spin 1/2 nucleus 
aiso measured in its rest- frame. Similarly, denotes the muon normalized polarization vector 
in its rest-frame, while A = ±1 is the massless Dirac^ neutrino helicity. 

Finally, the capture distribution is given by the expression. 



\Mx\\ (2) 



Here, dVtp is the element of solid angle associated to the outgoing nucleus of momentum p, Ai\ 
is the capture amplitude associated to a neutrino of helicity A, tpdO) is the IS state muonic atom 
Coulomb wave function measured at the origin, and C is a reduction factor which accounts for 
the effects of the nuclei finite size through the overlap of the muon and neutrino wave functions 
with the different nuclear electric and electroweak charge distributions of finite spatial extent. 
The detailed evaluation of this reduction factor C is discussed in Appendix 1 both for ^He and 
for the proton. 



2.2 The amplitude 

The parametrization for the effective interaction associated to muon capture, at the level of the 
u and d quarks, is taken to be, 



8M2 '^"'^ 



E 



+ 



T 

Vim 



VnO^'^PrnlJLdUayP^ 



~V2 



u 



(3) 



Here, P± = (l±75)/2 are the chirality projectors, and g, M and Vud are arbitrary real parameters 
which in the limit of the SM reduce to those of that Model, namely M = Mw, ^^/SM^ = 
Gp/V^ and Vud = cos 9c being the Cabibbo-Kobayashi-Maskawa (CKM) quark fiavour mixing 
matrix element, in which case all coefficients also vanish except for = 1. Finally, 

the coefficients h^l^^ ^"^^ arbitrary complex coefficients associated to vector, scalar and tensor 
interactions, with ryi (resp. 772) being the muon (resp. d quark) chirality, equal to the neutrino 

*In practice, there is no difference between a massless Dirac or Majorana neutrino in tfie case of purely 
{y — A) or (V + A) interactions, but the distinction becomes relevant as soon as other interactions are turned on, 
as done in this paper. By this assumption of a massless Dirac neutrino, we exclude the possibility of interference 
contributions between processes in which either the neutrino spinor field v{x) or its charge conjugate v'^{x) would 
couple in the amplitudes to the same quarks and leptons rather than their antiparticles, namely it is assumed 
that the muon leptonic flavour is conserved in the massless neutrino limit. 
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(resp. u quark) one for vector interactions and opposite to it for scalar and tensor interactions 
(the "*" symbol denotes complex conjugation throughout). Finally, without loss of generality, 
one may set h'^_^_ = = h^_, because of the identity = ie^'^^'^ (jptj /2 which implies that 

the terms multiplied by /i!^^ and simply vanish identically. 

This choice of parametrization is inspired by the one used (in the charge exchange form) 

S V T 

for muon decay in terms of coefficients g^j^' , with the first (resp. second) lower index being the 

chirality of the electron (resp. muon) 1331] . Note that a similar effective four-fermi interaction 

— S V T 

may be given for /3-decay in terms of coefficients f±±' with the electron then playing the role 
of the muon in (^) . 

To express the amplitude for nuclear muon capture in terms of the above parametrization, 
one also requires the hadronic matrix elements of the relevant quark operators in terms of the 
nuclear bound states. For this purpose in the case of a spin 1/2 isodoublet, let us introduce the 
spinor ipi for the Dirac field of the spin 1/2 nucleus on which the muon is captured, while 'tp2 is 
the Dirac spinor for the final spin 1/2 nucleus. In momentum space, with = P2 — Pi being 
the momentum transfer of the process and P2 (resp. p'^) being the momentum of the final (resp. 
initial) nucleus, we have the following parametrization in terms of g^-dependent form factors. 



< 2\d^^u\l >= V2 [^Fv-f^ + iFMa^u^ + Fsj^j Vi, (4) 

< 2|d7/,75u|l >= tlj2 ^-^47^75 + ^P^5^ + ^^'^^^'''^^2M^ 

<2\du\l>=W2{Gs)^i, (6) 

< 2|d75u|l >= V^(Gp75)'0i, (7) 

< 2\da^yu\l >= ip2 {GrfJ^u) ipi , < 2\da fj,uTou\l >= "02 (Gtct^,. 75) V'l, (8) 



where |1 > and |2 > denote the initial and final quantum nuclear states, respectively. 

This notation for form factors is the one used for the ^He-^H case, while it is more con- 
ventional in the proton-neutron case to denote the form factors for the vector and axial quark 
currents as gy, givi-, Os, dA-, Op and gr, respectively for Fy, Fm, Fs, Fa, Fp and Ft- We shall 
thus follow that convention in the case of the proton, and keep nevertheless the notation Gs, Gp 
and Gt for the form factors associated to the scalar, pseudoscalar and tensor quark operators 
even in the case of the proton. However, all expressions for observables in this section will be 
given in terms of the Fv^m,s,a,p,t form factors. Note also that the contribution of the induced 
pseudoscalar form factor Fp has conventionally been normalized to the muon mass, while those 
of all other recoil order form factors are normalized with respect to twice the mean nuclear mass 
2M. 

In the case of T-invariant interactions, both the effective coupling coefficients h^^' as 
well as the nuclear form factors are all real quantities under complex conjugation. However in 
the present analysis, this restriction is not imposed, and all these parameters are assumed to 
take a priori complex values. 

For the vector and axial current quark operators, one distinguishes first- and second-class 
form factors, Fy, Fm, Fa and Fp in the first case, and Fs, Ft in the second case. On the 
basis of CVC, the values for the electroweak vector and induced magnetic form factors may 
be related to those of the electromagnetic electric and magnetic form factors of ^He and 3H, 
or of the proton and the neutron. The value of the axial form factor Fa at zero momentum 
transfer follows from the /3-decay rate either of ^H or of the neutron. Its value at non vanishing 
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momenta transfers requires knowledge of its g^-dependency inferred from neutrino or pion elec- 
troproduction experiments. Finally, through PCAC, the value of the last first-class form factor, 
the induced pseudoscalar one, may be expressed as0 

j7PCAC/^2n _ 

^'i> = 7:;2 — ^2 — ;;2 — ^2—' y^) 



where m-,^ is the vr^ mass, f-n its decay constant, and g-nNio^) its nuclear coupling to the two 
nuclear states of masses Mi and M2. In the limit of exact isospin symmetry, which implies 
also exact G-parity invariance, the values for the second-class form factors Fg and Ft may be 
shown to vanish identically. Hence, one expects deviations from zero for these form factors 
(normalized to Fy or Fa) of only a few percent, as given by the ratio (m^ — ?7i„)/Aqcd of u and 
d quark masses to the QCD scale for example, or the ratio (m„ — rup) jirrin -|- m^) in terms of the 
neutron and proton masses. Bag model or QCD sum rule evaluations of these form factors in 



the case of the (proton,neutron) doublet do indeed bear out such an expectation |P6|, 37, 38, 39 1. 
In the same manner, one could wonder about the effects of isospin symmetry breaking for 
the Fy and Fm form factors obtained through CVC. In this particular case, thanks to the 
Ademollo-Gatto theorempO|, isospin breaking corrections are only of quadratic order in the 
ratio (mrf — m„)/AQCD) hence negligeable for our purposes. Let us also recall that CVC has 
been very well established through precision studies in /?-decay|4l|], and that even though limits 
exist on second-class form factors from correlation experiments in /3-decay, the stringency of 



these constraints are nowhere close to the theoretically expected values for gs or 5(t[33|. 

In the matrix elements for the scalar, pseudoscalar and tensor quark operators, we have 
chosen not to include recoil order induced contributions, for the reasons mentioned already in 
the Introduction. The values for the genuine scalar, pseudoscalar and tensor form factors G5, 
Gp and Gt, cannot be inferred from any experiment yet. One thus has to rely on specific 
model calculations for QCD dynamics, such as the bag model or QCD sum rules. However, no 
such results are available at present, and one may only reasonably guess that these form factors 
should take values on the order of unity, within a factor which at worst could be of order ten. 

Hence, even though the parametrization of the nuclear matrix elements of the relevant 
quark operators in terms of form factors is only a phenomenological representation of our ig- 
norance of the microscopic nonperturbative quark dynamics, this approach to nuclear muon 
capture allows nevertheless for explicit predictions independently of the details of nuclear mo- 
dels, relying only on the results of other experiments and the power of symmetry principles [Q]. 

In terms of this parametrization of the nuclear state matrix elements, the effective muon 
capture amplitude at the nuclear level is given by, 



,2 



8M2 



(/i^i,,2)*i^/.7^(l + f?i75)MX 



X '^2\lii,{Fy r]2FAlh) + -^{Fs + ri2—Fpj5) + i<y,j,u^{FM + V2FT-f5Ml + 

2M m,, 2M 



+ « V2 ) + ^175)^ i'2iGs - mGp^5)i'i + 



(10) 



To be precise, these expressions for Fp assume implicitly that both q -dependencies for -Fa(? ) and gnN{q ) 
are identical Q]. Note also that through relations such as these, any experiment leading to a precise value for Fp 
would also imply a precise value for the associated pion- nucleus coupling constant |B5|. 
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2.3 The method of calculation 



The remainder of the calculation requires now the evaluation of |A^aP- If one were to proceed by 
"brute force" , using the usual trace techniques for such calculations, one would quickly run into 
unmanageable expressions, because of the many contributions stemming from all the interactions 
represented in the amplitude A^a- Actually, it is possible in the present case to take advantage 
of the fact that the initial system is at rest, and that to a very good approximation the initial 
muon and nucleus may be considered to be also at rest]]. 

For that purpose, it turns out that the Dirac representation of the Dirac-Clifford algebra of 
7^^ matrices is the best suited for the problem. Plane wave solutions to the massive Dirac equa- 
tion are then of the following form, for "positive" and "negative" energy solutions, respectively 
(our Minkowski metric signature convention is (H )), 

(11) 

where x±(s) are bi-spinors such that 

^ ■ s X±{s) = ± X±{s) , X±{s)x±{s) = I, (12) 

with fj* (i = 1, 2, 3) the usual Pauli matrices, while s is a unit vector in three dimensions which 
in fact corresponds to the spin polarization vector of the particle (see Appendix 2). Further 
properties of these bi-spinors are 

X+{-s) = ix-{s) , x^{-s) =ix+{s), (13) 

Xv{s)xl{s) = ^[l+w-s] , r? = ±. (14) 

In fact, it is only the latter relation which is of essential use in the calculation of |A^aP- 

Even though this is not important for that calculation, it is also possible to give the explicit 
expressions for these bi-spinors. Associated to the spherical angular parametrization of the unit 
vector s. 



s 



I 



\ 



sm f cos ip 

sin0sin(/? I , (15) 



cos ( 



one has 



... / e-^^/2cos0/2 \ / -e-*'^/2sine/2 \ 

^+^')=l e^^/2sin0/2 ) ' ^-(«)=l e^./2cos0/2 j" ^'^^ 



Finally in the case of solutions to the massless Dirac equation, "positive" and "negative" 
energy plane wave solutions of helicity A = ±1 and momentum k are given by 

where of course k = k/\k\. 

^This assumption amounts to ignoring those small relativistic corrections which are related to the velocities of 
the bound muon and nucleus, which in the present case is indeed totally justified. A fully satisfactory relativistic 
treatment of such a bound state problem in QED is still not available, and cannot be used here to estimate small 
corrections which in any case will be at most of order (q^)^, namely the squared velocity of the bound muon. 
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In order to apply these considerations to the capture amphtude (|lCl|) , let us denote by 
X^^is^J.)^ Xi(si), X2{s2) and Xx{~p) the bi-spinors associated to the muon, the initial and final 
nuclei, and the neutrino, respectively (for the first three states, they thus correspond to bi- 
spinors of x+{s) type). After substitution in (|lO|), the capture amplitude associated to a massles 
neutrino of helicity A then reduces to the expression, 
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8M2 



Vud V^mMi. 



2^s 



X 2< x{Xf, xl 



i^f + (<)-<))p.a 



Xl+ 



+ A X^i xl 



Hl'^'iip X ay - Hp a' 



Xl 



with coefficients H 



(5,p,y,A) 

A 



defined by the following relations 



(18) 



H 



H 



(S) 



J2 "l^f^xm 



H 



(V) 



+ I E ^2 A,. 



Gp^ - 2X{htxx) G)j: 



I r?2=± 



(19) 

(20) 
(21) 



E ^2/lA^2 



.(2) 



in which the following combinations of form factors are introduced. 



G^v^ = {V~s-M2) 



^ 2M *^ 



+ (^/s + M2) 



Fv + ^-^^F, 
^ 2M * 



(22) 



(23) 



G 



(1) 



s-Mo 



G 



(2) 



G 



(2) 



(V^-M2) 

: (Vi + M2) 



Fy + 



+ (Vs + M2 

Ml +M2 



2M 



2M 



Fm 



G 



(1) 



Ml - M2 ■ 

F4 -Ft 

2M . 

i+M2)Gs, 



G« = (Vi-M2)Gp, 



G 

G, 



(1) 

T 
(2) 



S - M2)Gt, 
S + M2)Gt. 



(24) 

(25) 

(26) 

(27) 
(28) 
(29) 
(30) 



The calculation of |A^aP then proceeds using the trace properties in (|T^. The advantage 
of the above approach is that the combinations of form factors and coupling coefficients which 
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are relevant appear from the start, ever before proceeding to the calculation of traces. Were one 
to first calculate the traces as is usual, the task would quickly become impossible in the general 
case considered here (presumably, this is the reason why only numerical expressions were given in 
Ref . , even though the situation considered there was far less general) . Still, using the present 
approach, the length of the calculation is of some importance. Note also that no nonrelativistic 
expansion in the amplitude is effected at any stage of the calculation, in contradistinction to all 
other analyses which are based on Ref.Q] in which such an expansion in 1/M is indeed applied. 

Finally, the capture distribution is thus given by, 

dV^ _ |V(0)|2 / 52 \2 ^2 



d9-p 327r2 l8M2 



Kd-^A, (31) 



with, 

Rx = 

+ (l-Xp-Si) {Cf) (l + (p.5^)(p.s2))+Z)J) ((|^.52)-(p.S^)(p.S2))+4'^p.(s^XS2)| + 

+ (l + Ap.S2) {Cf (l-(p-s^)(p-Ji)) + Z)f (^^s^.s,)-{p.s^){p-si)) + E^^^p-is^x.h)^, 

(32) 

in which the following final definitions apply, 

cM ^ |^(5)|2 ^ |^(P)|2 _ 2 |^(^) + A/7(^)|2, (33) 

= -2ARe (i/f (<V) + 2 \H^P + AfJ^^I^, (34) 
= +2ARe (Ff (<V) + 2 \H^P + A<)|^ (35) 

^(A«)^|^{S)|2_|^(P)|2^ (36) 

I?!') = -2Re ((/^f - A//f ))(/7f ) + A^i^))*) , (37) 
^ = -2 Re ((//f ^ + XhP){H^P + Af|^^)*) , (38) 
=2Im(//f (FfV), (39) 



= -2AIm ((/^f - A<))(<) + A<))*) , (40) 
E^^^ = -2AIm ((//f ^ + XhP){H^P + A-ffJ^V) • (41) 

Note that these results show that T-odd effects can only appear through the triple cor- 
relation coefficients of E^^'^'"^^ type, related to contributions which are pure imaginary under 
complex conjugation and involving necessarily always at least two polarization vectors. 
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2.4 Capture distributions and final state polarization 



The calculation thus leads to the capture distribution as given in (|32|). For a massless neutrino 
of given helicity A {i.e. a Weyl neutrino, or a neutrino whose helicity is measured!), (|32|) gives 
the final expression. For a massless Dirac neutrino produced with either helicity (which is then 
not measured), one needs still to sum the result over A = ±1. Finally, and independently of the 
neutrino sector, in order to obtain the unpolarised nuclear final state distribution, one simply 
needs to sum the two results obtained from (^) for §2 = sq and §2 = —sq. Since the result ( p2[ ) 
is at most linear in §2, this amounts to setting S2 = and to multiply the result (^) by a factor 
two. 

In general, (|3^ ) thus has the parametrization, 

R^=Mx + S2- -Dx, (42) 

where the expressions for the coefficients J\f\ and may be read off directly from (^). The 
capture distribution of the final state nucleus is thus given by, 



dTx 
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8M2 



(43) 



while the final state spin 1/2 nucleus has the following average polarization vector, 

1 ^ 

< S2 >x= TT ^A- 

A/a 



(44) 



For example taking an initially unpolarized state with 
tion of the final nucleus is given by 



= si, the average polariza- 



(unpolarized) 
< «2 >\ 



A- 



C 



(2) 



(1) 



.(2) 



(45) 



Hence, a transversally polarized final nucleus requires an initially polarized muonic atom. 

However, all the results discussed so far refer to polarization states defined in terms of 
the individual polarization vectors s^, si and §2, rather than the hyperfine states {8,171) of the 
initial muonic atom, with (S" = l,m = ±1,0) and {S = 0,m = 0). As explained in Appendix 2, 
it is of course possible to use the above results to determine the final state distributions and 
polarizations associated to each of these hyperfine states. Since it is difficult to imagine how the 
final state polarization of the neutron or ^H could ever be measured to any degree of precision, 
here only the relevant expressions for the hyperfine capture distributions are presented. 

Given the hyperfine states 5 = 0,1 and their projections m = and m = 0, ±1 along 
some arbitrary quantization axis, the associated hyperfine capture distributions are given by. 



dr 



(S,m) 



\m\' 

167r2 I 8M2 
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(S,m) 



(46) 



with 



R 



(i,±i) 



(1) 



R 



(1,0) 



(2) 



+ D 



(2)> 



^X{c[^Uc\'')P^{co^ 



D 



{2)^ 



(Cf>+D 



PofcOSi 



{2)' 



P2{cose), 

(47) 
(48) 
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R 



(0,0) 



(Ci'nc«) + 2(cf 



D 



(2)n 



AD 



where the usual Legendre polynomials are 

Pi (cos 6) = cos 9 



P2{cos6) = - cos^ 



(2) 
A ' 

1 

2' 



(49) 



(50) 



while 9 is the angle between the spin quantization axis and the normalized momentum vector p 
of the outgoing spin 1/2 nucleus. 

In particular, the statistical capture distribution is given as in (^6[) with the coefficient 



^(S',m) ^j^gj^ obtained from 



c(^) + C['' + c 



(1) 



.(2) 



\H 



R 



(5) 1 2 



(1,-1) 



+ + R 



(0,0) 



{^')|2 



+ 2|F 



XH 



(^)|2 



(51) 



a results which is indeed ^-independent as it should. 

Integrating these distributions leads to the associated capture rates, A^, Xj^ and Aj'^* for 
the singlet (S), triplet (T) and statistical rates, respectively, of the form 



A 



S,T,stat 
A 



i-TT 



9 



8M2 



V, 



ud~ 



R 



S,T.stat 
A ' 



(52) 



with i?!*'^* for the statistical capture rate given in (^Tj) already, while for the singlet and triplet 
capture rates, one has, respectively. 



c 



(2) 



R 



(m) 

A 



2D 

(2) 



(2) 



D 



R 

(2) 



(0,0) 
A ' 



It is also possible to represent the hyperfine capture distribution (|4 
formg, I 



1 , (stat) 



1 + A^Pa + A^P^ cos e + A^Pti^ cos^ ( 



(53) 
(54) 

in the following 
(55) 



where A'^, A^ and A^ are specific coefficients, the latter two known as the vector and tensor 
analyzing powers of the final state nucleus, respectively, while A^ is a measure of the hyperfine 
effect on the statistical capture rate since one has 



A! 



-^A 



A! 



4 AJ^* 



3X1 +xr 



(56) 



Finally in (|55D, the coefficients PA,v,t are the following combinations of the hyperfine populations 

Ns,m 

Pa = iVi,i + Nifl + Afi,_i - 3iVo,o = 1 - 4iVo,o, (57) 
= iVi,i - Afi,_i , Pt = iVi,i + iVi,_i - 2A^i,o, (58) 
such that Ni^ + A^i _i + A^i^o + A^o.o = 1- 

In terms of the quantities introduced above, one then finds. 



^A 



IC 



(2) 



2D 



(2) 



^stat 



A. 



X{C[^'> + C 



(l)^ 



2 + Zjf 

2 ^stat 
A 



(59) 



Note that when observables are considered in which the summation over the two neutrino helicity 
states A = ±1 has been effected, this summation has to be applied separately in the numerator 
and the denominator of each of the expressions given in this subsection. 
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3 The ^He Case 



3.1 Physical inputs 

The basic kinematical input used in the case of muon capture on ^He is as fohows, 

Ml = 2808.392 MeV , Ma = 2808.928 MeV , ^/s = 2914.039 MeV, 
u = 103.22 MeV , uo = 2810.82 MeV , to - M2 = 1.90 MeV. 



(60) 



As explained in Appendix 1, the overlap reduction factor C in this case takes the value 
C = 0.979|2|, ^, while the latest value quoted in Ref.|33] for the CKM ud mixing angle is used, 
namely^ 

Vud = 0.9750 ± 0.0008. (61) 

This value results from a global fit which includes the constraints of unitarity of the CKM mixing 
matrix, rather than the value which would follow solely from the 0~*'-0~'~ superallowed /3-decay 
/t- values. 

Values for the nuclear form factors at the relevant momentum transfer qf = — 0.954m^ 
must also be specified. For this purpose, we refer to the discussion in Refs.[§, ^ and use the 
values advocated by these authors. For the electroweak vector and induced magnetic form 
factors, one has 

Fviqf) = 0.834 ± 0.011 , FM{ql) = -13.969 ± 0.052. (62) 

The situation for F^iqi) is somewhat more delicate, since only its value at = is known from 
/9-decay of '^H, while its g^-dependency may only be inferred through a nuclear model calculation 
of the associated mean square charge radius [§, ^. This leads to the following value 

FA(g?) = -1.052 ± [0.005 - 0.01], (63) 

where the interval for the uncertainty is an attempt to reflect the lack of precise knowledge of this 
form factor. The lower uncertainty of 0.005 stems from the uncertainties on the /3-decay rate, 
while the corrections due to mesonic exchange currents may only be estimated in the nuclear 
model calculation, leading to a total uncertainty of 0.007 on FA(gi)[^, Thus the bracket 
[0.005 — 0.01] represents a conservative evaluation of the uncertainty on that form factor, which 
will be carried along throughout our analysis later on. As it turns out, the statistical capture rate 
is rather sensitive to that quantity, so that any improvement on the evaluation of its uncertainty 
leads to an improvement on the values of other quantities that one infers from experiment. 
However, it is difficult to see how such an improvement on F^iqi) could be achieved in practice, 
since it would require a measurement of the g^-dependency of the ^He-^H axial form factor 
FA{q^). 

When it has to be specified, the value used for the induced pseudoscalar form factor Fp 
is that inferred from the PC AC relation (^), which corresponds to 

Fl^^^{ql) = -20.73 ± [0.10 - 0.20], (64) 

where the uncertainty range indicated in brackets thus corresponds to the associated range in 
the value used for F^(g^). 

*The uncertainty on this value is irrelevant for the '^He case , but must be included in the prediction for muon 



capture on the proton at the precision level reached in Sect. 4.2. Note also that this value for V^d differs somewhat 
from that used previously |l^, |l^], which implies that the numbers quoted here differ somewhat from those 
quoted earlier, but in no physically significant way whatsoever. 
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Finally, the remaining two second-class form factors for the vector and axial quark ope- 
rators, Fs{qi) and FT^qi)-, are taken to be vanishing, even though for the reasons explained in 
Sect j2 .21 , one should expect that their absolute values, normalized to those for Fv{qi) and ^4(5^), 
respectively, would be on the order of 0.02. As will become clear later on, the statistical capture 
rate is in fact rather insensitive to these two form factors, so that this reasonable approximation 
is certainly sufficient. 

These are all the nuclear form factors required when considering muon capture within 

the SM. Beyond that model however, the genuine nuclear scalar, pseudoscalar and tensor form 

factors, Gsiqf), Gp{qD and Griqi), ^^^o required. As mentioned previously, values for these 

quantities are not known, but are expected to be on the order of unity within a factor of at most 

S V T 

ten. Nevertheless, any constraints put on the effective coupling coefficients h^^' beyond the 
SM will thus involve at present these nuclear form factors as well. 



3.2 Within the Standard Model 

Capture within the SM corresponds to the amplitude Aix in (10) with only the = 1 effective 



coupling coefficient turned on. Given the above numerical values, as well as the general expres- 



sions of Sect. 2.4, it is straightforward to determine theoretical values for all the quantities which 
enter the general final state triton distribution associated to hyperfine states as parametrized in 
(H). Onefind^g, 

= 1929 ± [28 - 46] s-\ 
A"^ = 1351 ± [16 - 19] s-i, 
^stat ^ 1496 ± [12 - 21] s-\ , . 

A A = -0.0967 ± [0.0061 - 0.0069], ^ ' 

A = +0.524 ± [0.0057 - 0.0061], 
At = -0.3790 ± [0.00074 - 0.00123], 

where each time the indicated uncertainty includes all uncertainties of all the theoretical input, 
while the range indicated by the bracket corresponds to the range implied by the uncertainty 
on ^^(^1). Obviously, some of these results coincide with the specific quantities also computed 
in Refs.]Q, ^ in this particular case. 

The disparity in the uncertainty ranges for these quantities stems from their sensitivity to 
all theoretical inputs. This sensitivity of an observable O to a parameter Fx may be specified 
in terms of the variations, evaluated in the SM, 

a{0;Fx) = ^^ if Fx ^ ; a{0;Fx) = ^^ if Fx = . (66) 

The values of these quantities for all the above hyperfine observables are shown in Table |l[ Note 
that except for A^, all these observables are rather sensitive to the axial form factor Fyi(g^), 
and thus to its bracket of uncertainty. Sensitivities to Fp^qf) include the results of Refs.]|2|, |^, 
showing in fact that except for the triplet capture rate, the statistical rate which has been 
measured [^] is the least sensitive observable to that form factor, with the obvious drawback 
that all other observables are still far more difficult to measure. 

Given the experimental result for the statistical capture rate|12]. 



Xf^^ = U96 ± 4 s-\ (67) 



^Note that the hehcity index A is suppressed, since whether one sums or not over the hehcities of the final 
neutrino which is only of left-handed chirality is irrelevant in the massless limit in the present case. 
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it is possible to infer a value for any given parameter once the values for all other quantities are 



specified. Thus given the inputs discussed in Sect.3J. except for the value for Fp(g^), the result 
® implies|l2|, 0, |15|, 

Fp{qf) = -20.69 ± [1.57 - 2.74][exp : 0.48], (68) 

where the first uncertainty bracket includes both the experimental error as well as all the un- 
certainties on the input form factors, the range corresponding again to the range in the FaCqi) 
uncertainty, while the last number in brackets represents the uncertainty following only from 
the experimental error in (|67|), namely without any of the errors on the other inputs. This 
last number thus indicates the range of improvement that could be achieved by reducing the 
uncertainties on the input form factors Fv{qi), i^Af (^i) and Fa (91)1 and especially on the latter 
one. Compared to the PCAC expected value for this induced pseudoscalar form factor, the 
above result is thus in confirmation of the PCAC prediction with a precision ranging from 8% to 
13%. In order to translate this conclusion in terms of the nucleon form factor gp, it is necessary 
to include in the nuclear model calculation all meson exchange corrections and their associated 
uncertainties, thereby leading to a 19% precise test of PCAC with a value in agreement with the 
prediction |T^. Note also that given the PCAC relation (^), the above value for Fp(gf) inferred 
from experiment may in turn be used to determine the tt^— ^He-^H nuclear coupling constant 
to much improved precision j3 5(1 . 

A combined fit to two independent observables measured to great precision, such as A^*'^* 
and the vector analyzing power Ay, would allow a model independent determination of both 
FA{qi) and Fp(g^). However, it may be shown that in order to obtain a result for FA{q\) with 
an uncertainty less than the present value of 0.01, would require a measurement of to better 
than 1%, no small feat indeed! In particular, when used on its own, a 1% precise measurement 
of Ay, centered onto its theoretical prediction, would imply the following uncertainty range for 



Fp{q1) when using for the other form factors the values quoted in Sect. 3.1, 

Fp{ql) : [0.80 - 0.80] [exp : 0.55]. (69) 

Turning now to the second-class form factors Fs^ql) and Fxiql), using the PCAC 
prediction (64) for Fp{q\), the experimental result (B^) implies either 



Fs{qi) = 0.026 ± [1.17 - 2.02] [exp : 0.38], (70) 

or 

Friql) = -0.031 ± [1.42 - 2.45] [exp : 0.46], (71) 

when either form factor is turned on while the other is still set equal to zero. Due to the small 
sensitivity of the statistical capture rate to these two parameters, these constraints are thus 
extremely poor, but nevertheless they improve somewhat the situation existing in terms of the 



nucleon form factors gs and 5t[23] if one is willing to extrapolate without correction from the 
hydrogen case. This conclusion would not be much improved were a 1% precise measurement 
of Ay to become available, since the corresponding uncertainties ranges are then, 

Fs{ql) : ±[0.85 - 0.89] [exp : 0.58] ; Friql) ■ ±[0.80 - 0.83] [exp : 0.54]. (72) 

Similar considerations could of course be developed on the basis of other observables still, but 
we shall refrain from doing so since they are beyond the reach of experiment at present. 

Finally, let us point out that if Fp{qf) is allowed to vary within its uncertainty bracket in 
(|68|), the corresponding values for Fs{qi) and F'j'{qi) then also vary accordingly, but still within 
their respective uncertainty brackets. Hence any dependency of the results for Fs^t on the value 
assumed for Fp is consistent with their own present uncertainties. 
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3.3 Beyond the Standard Model 



Let us now turn to the reach for physics beyond the SM offered by the experimental result 
(pTl). Defining sensitivities a{0;hx) of observables to the effective coupling coefficients h^^' 
in the same manner as in (^) with respect to the vanishing second-class form factors Fs^t, the 
corresponding results are given in Table 0, using as input for nuclear form factors the values 



discussed in Sect.^3| as well as unit values for the scalar, pseudoscalar and tensor form factors 
Gs,p,T- In fact, these observables are not sensitive (in linear order) to the coefficients h^^., 
and because of the fact that only a left-handed neutrino couples to muon capture in the 
limit of the SM. Note that some of these observables are quite sensitive to the tensor coupling 







hi 






hi 


hi 


^hl_ 


+ hl 


hi 


= hi_ 


-hi 



/iJ^_/2, for which stringent constraints will thus be inferred from (67). 

Indeed, the experimental statistical capture rate ( [67| ) implies the following results, turning 
on only one coupling at a time, 

= 0.00049 ± [0.0224 - 0.0385] [exp : 0.00723], 

^ = 0.00047 ±[0.0218 - 0.0374] [exp: 0.00702], 

^ = 0.00048 ± [0.0221 - 0.0379] [exp : 0.00712], 

^ = -0.0322 ± [1.49 - 2.55] [exp : 0.48], (73) 

\h%_ = 0.000031 ± [0.00143 - 0.00245] [exp : 0.00046], 

h^_ = 1.00009 ± [0.00415 - 0.00712] [exp : 0.00134], 

h'^^ = -0.000222 ± [0.0102 - 0.0176] [exp : 0.0033]. 

These values assume implicitly that the form factors Gs,p,t have been set equal to unity. Note 
also that the coefficients hi and h^ do define actual scalar and pseudoscalar interactions at the 
quark-lepton level, for a neutrino of left-handed chirality. 

Again, it may be checked that allowing Fp to vary within its uncertainty bracket (|68|), 
each of the above coefficients then also varies essentially within its own uncertainty bracket. 
As could be expected from Table the results for the scalar couplings hl_, hl^ and hi are 
already quite stringent, and in fact improve present limits on such couplings both in the muonic 
as well as in electronic sectors. However, the most satisfactory result is undoubtedly obtained 
for the tensor coupling h'^_, which is brought down into the per mille level. Of course, definite 
conclusions as to the actual limits for such physics beyond the SM implied by the experimental 
result (^) would require the evaluation of the form factors Gs{qi), Gp{qi) and Gt(?i)- The 
above limits on are translated in what may be physically more meaningful terms in Sect]^. 

Again for comparison, it is interesting to establish what a 1% precise measurement of 
Ay, centered onto its theoretical prediction, would imply for the same set of effective coupling 
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coefficients. Correspondingly, one finds for the associated uncertainty brackets: 



hl_ 


±[0.0164 


- 0.0170] [exp : 


0.0111] 




±[0.0157 


- 0.0164] [exp : 


0.0107] 


hi 


±[0.0161 


- 0.0167] [exp : 


0.0109] 


K 


±[0.81 - 


0.84] [exp : 0.55], 




±[0.0201 


- 0.0210] [exp : 


0.0136] 




±[0.0144 


- 0.0150] [exp : 


0.0098] 



(74) 



Hence, such results would improve somewhat on the situation for some coefficients in ( [73| 
following from the experimental statistical capture rate 



4 The Hydrogen Case 
4.1 Physical inputs 

In the case of muon capture on the proton, the kinematic variables are such that we have 

= 1043.93 MeV , i/ = 99.146 MeV , w = 944.78 Me V , a; - Ma = 5.22 MeV. (75) 

We also use for the CKM ud matrix element the value for V^d quoted in (|6T1), including its 
uncertainty. Indeed, the latter contributes 0.16% to the final uncertainty on the theoretical 
prediction, and should thus be included given the aim of a 0.5% precise measurement of the 
singlet capture rate[|l^, |^]. As to the overlap reduction factor C, the value used is that detailed 
in Appendix 1 in the case of the proton, namely C = 0.9956. 

Let us turn to the issue of the nucleon electroweak form factors, to be evaluated at the 
momentum transfer = — 0.877m^, beginning with the vector ones, gviQo) and gAiiQo)- Since 
through CVC, the values for these nuclear form factors are related to the electromagnetic ones 
for the proton and the neutron whose charge radii are very well known[42|, it becomes possible 



to infer very precise values for gvioi) and guilo)- From the values rf = 0.765(1 ± 0.01) fm and 
r2 = 0.893(1 ± 0.01) fm (see Ref. [p2| for the meaning of these parameters), as well as the proton 
and neutron anomalous magnetic moments and the fact that gv{Q^ = 0) = 1, the following 
values apply 

gy{ql) = 0.9755 ± 0.0005 , gniql) = 3.5821 ± 0.0025. (76) 

The value for the axial form factor 5a (^o) inferred from the neutron decay rate, which 
implies |p3|, 

gj^{q^ = 0) = 1.2670 ± 0.0035, (77) 
as well as the mean square axial charge radius [|^, § |9| 

rA = 0.65 ± 0.03 fm. (78) 

Consequently, one finds 

gA{ql) = 1.245 ± 0.004. (79) 
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Note that since the capture rates are rather sensitive to the axial form factor (see below) , if the 
value for gA{(f = 0) were to change again, the values for capture rates would have to be adapted 
appropriately. 

The value for the induced pseudoscalar form factor gp{qQ) is given by the prediction based 
on the chiral symmetries of QCD||, H, 

9pi%) = ^2 2 o9A[0)mf,Mrj^, (80) 

where f-,^ = 92.5 ±0.2 MeV is the usual pion decay constant and g-j^NN the pion-nucleon coupling 
constant. The latter quantity has been the issue of much debate, but using the recent precise 



estimate of Ref . [ 44 1 , 



y^^^ = 13.37 ±0.09, (81) 

one finds 

c,^(g2) = 8.475 -t 0.076, (82) 

namely a chiral symmetry prediction with a precision better than 1%. 

The second-class form factors gs and gx are taken to be identically vanishing, an approxi- 
mation discussed in Sect.^]^ which is justified in the limit of exact isospin symmetry. However, 



the effects of isospin breaking are estimated to be small|36, 37, 38, 39], on the order of 0.02 for 
\gs/gv\ and Igx/gAl, and given the sensitivity of the capture rates to these two form factors (see 
below), their contribution may safely be neglected. 

Finally, there is the issue of the nucleon scalar, pseudoscalar and tensor form factors 
Gs{qo), Gp{qQ) and Gxiqa), whose values are unknown at present. Any estimate would of 
course be welcome, but the expected result should be of the order of unity, within may be a 
factor of at most ten, as mentioned previously. 

4.2 Within the Standard Model 



Given the different inputs discussed in Sect. 4.1, it is straightforward to consider predictions for 



the observables in muon capture on the proton within the SM. Setting all effective coupling 
coefficients to zero except for = 1, one then finds, 

AS = 688.4 ± 3.8 s-\ 
A"^ = 12.01 ± 0.12 s-\ 

;^stat ^ 181.11 ±0.98 S-\ , . 

Aa = -0.9337 ± 0.0007, ^ ^ 

= 0.00371 ± 0.00030, 
At = -0.06260 ± 0.00052. 

Note the small values for the latter two observables, rendering their experimental determination 
essentially impossible, given the small polarizations available for the muonic hydrogen atom. 

As in the case of muon capture on ^He, it is interesting to consider the sensitivity of each 
of these observables to the different input nucleon form factors, in particular gpiq^). These 
sensitivities are defined as in (^6|), and their values are presented in Table ^. Note that with 
respect to the measured capture rates, namely the statistical one for ^He and the singlet one for 
the proton, their sensitivity to the induced pseudoscalar form factor is essentially comparable, 
as is also the case for the other nucleon form factors. 
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In fact, assuming the goal reached |TJ] of measuring the singlet rate on the proton to a 
precision of 0.5% with a central value equal to the above theoretical prediction precise to 0.55%, 
one may infer the following value for the induced pseudoscalar form factor, 

9p{<ll) = 9p{ql) ± 0.327 [exp : 0.230], (84) 

where the second indicated uncertainty follows only from the experimental precision of 0.5%, 
while the first also includes all the errors on the theoretical inputs for the other form factors 
and the CKM matrix element Vud- In other words, a 0.5% precise measurement of the singlet 
capture rate implies a 3.9% precise determination of and for the test of the chiral symmetry 
prediction of that value. In the same way as in Ref . [p5|] , the same result may also be used to infer 
a value for the pion-nucleon coupling constant q-kNN: assuming that gp{qQ) takes its expected 
value (H), 

^^iVTv = 13.37 ±0.49, (85) 

hence a 3.7% precise determination of that quantity. 

Pursuing along the same lines, assuming the value for gpiq^) set by the chiral symmetry 
prediction, and a measurement of the singlet rate precise to 0.5% in the manner described above, 
each of the second-class form factors may also be constrained to the following precision, 

gs : ±0.314 [exp : 0.216] , gr : ±0.307 [exp : 0.210], (86) 

thus providing much of an improvement on the present situation [p3[], but still an order of mag- 
nitude away from the expected range of values. 



4.3 Beyond the Standard Model 

S VT 

The sensitivities of the considered observables to the effective coupling coefficients h^^' are 
given in Table ^ which displays some interesting differences with the '^He case, but again not 
between the singlet capture rate on the proton and the statistical rate for '^He. 

In order to assess the potential reach offered by the singlet rate for physics beyond the SM, 
let us again assume a 0.5% precise measurement of that observable centered onto its theoretical 
prediction. One then finds the following uncertainties, for each of the relevant effective coupling 
coefficients turning them on each one after the other. 







hl_ : 


±0.0168 [exp : 0.0115], 






hU : 


±0.0187 [exp : 0.0128], 


hi 


=hl_ 


_+hl^ : 


±0.0177 [exp : 0.0121], 


hi 


= ht 


- h^ 


±0.336 [exp : 0.230], 






2"-+- • 


±0.00140 [exp : 0.00096] 






h^_ : 


±0.00364 [exp : 0.00250] 






hV^ : 


±0.0113 [exp : 0.0065], 



in which it is implicitly assumed that the form factors Gs,p,t have all three been set to a unit 
value. Consequently, given the results in (|73D, the physics reach offered by the singlet muon 
capture rate on the proton measured to 0.5% precision improves by factors of at least two to 
three the results achieved already from the available statistical capture rate on ^He. 



18 



5 Specific Models Beyond the Standard Model 



The effective interaction (y) in terms of the couphng coefficients h^^' provides a model inde- 
pendent parametrization for any low-energy contributions stemming from whatever new physics 
there is lurking behind the confines of the Standard Model. However, taken as such, the con- 
straints (l?^ ) obtained from the experimental statistical capture on ^He probably do not mean 
much in terms of the energy scales or coupling strengths implied by such constraints. To develop 
an idea for the latter kind of data however, it becomes necessary to consider specific models for 
physics beyond the SM. Three general such classes will briefly be considered here, with the 
notations defined in Appendix 3. Other possibilities come to mind, such as for example models 
with large extra dimensions which may be amusing to investigate in the same vain. Also, when 
the singlet capture rate on the proton will have been measured, similar considerations may be 
developed as well, improving to some extent on the present results. 



5.1 Left-Right symmetric models 

Turning first to left-right symmetric models (LRSM) , only the limits on the couplings h^_ and 
hY^_^_ need to be considered, since only vector and axial interactions are implied in such models 
(up to small scalar Higgs exchange interactions which may be ignored, see Appendix 3). Given 
the value from (|7^), 

= 1.00009 ± 0.00712, (88) 

different considerations may be developed. First, one may view this result as a constraint on 
the universality of the electroweak interactions whose most stringent limit in the electron- 
muon flavour sector follows from vr^ decays|45, ^ at the 0.4% level. The above limit implies a 
universality constraint at the 1.4% level only (in terms of (/i^_)^). 

The above result for may also be viewed as a constraint on the unitarity properties 
of the CKM quark flavour mixing matrix, in terms of the ratio T^d/V^'*'''^^. Indeed, in muon 
capture one is dealing entirely with the muon flavour sector only, while V^d usually involves 
different /3-decay processes, hence the electron flavour sector. Therefore, another way to read 
the constraint on h^__ is to say that the experimental statistical capture rate on ^He confirms 
the unitary-constrained value for Vud, to a level better than what is achieved in terms of the 
0"'"-0"'" superallowed /3-decays where conclusions are somewhat dependent on the nuclear models 
used to evaluate the radiative corrections |Q. 

Finally, within LRSM, the above result for does not imply any limit on the mass 
of the extra charged gauge boson, since an expression similar to that given in Appendix 3 for 
the coefficient also applies to the similar coefficient fY- relevant to /3-decay, so that when 
reexpressing all couplings in terms of physical quantities (such as the muon decay rate, and so 
on), such factors essentially cancel. As a matter of fact, the sole contribution which survives this 
comparison of muon decay, /3-decay and muon capture amplitudes, is that in LRSM the CKM 
leptonic fiavour mixing matrices may be different for the electron and muon fiavours. Thus in 
fact, the above limit on hY_ translates into the following constraint on leptonic fiavour mixing 



in these two sectors, when the mixing angle ( of the charged gauge bosons is ignored, ( = 0|15|, 

{Vf, - Ve)vur'^S'^ = -0.00018 ± 0.01424, (89) 

with 
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in which the summation over the index i stands for ah neutrino mass eigenstates whose produc- 
tion is kinematically ahowed in muon capture on the one hand, and in /3-decay on the other (in 
which case, their mass is taken to be neghgeable as well). The other parameters are defined in 
Appendix 3. 

Similarly, the result on in ( [7^ ) may be translated in terms of parameters of LRSM 
as (see Appendix 3), 

rt Re (^e'^^Vud) = 0.000222 ± 0.0176. (91) 

In the case of manifest LRSM with r = 1, Vud = 1 and lo = 0, this result provides a limit on the 
mixing angle C, which does not improve such limits stemming already from /3-decay processes p3|. 

In fact, again in the limit of a vanishing mixing angle ("j the vector analyzing power yl-y 
reduces to 

^LRSM _ - r 5 Vf^Vu ,SM /Qr,^ 

lc=o- 1 + ' ^ ' 

with its SM value. However, even a measurement to 1% of Ay would not imply a lower 
bound on better than 260 GeV (95% C.L.) in the manifest LRSM. 

5.2 Contact interactions 

Let us now consider the possibility of contact interactions (see Appendix 3). Given the para- 



metrization of such interactions, it is a simple matter to translate the limits (73) in terms of 
the associated compositeness scales, with the following lower bounds (in an obvious notation) 
all given at the 95% C.L., 

Al_ > 1.60 TeV, 
Af+ > 1.62 TeV, 
A^ > 1.61 TeV, 

(93) 

A^ > 196 GeV, 
A^_ > 6.34 TeV, 
A^;^_ > 2.36 TeV. 

Clearly, the limits on A^_ and A!^_ are quite competitive with recent collider results in some 
of these channels p6|| . One also has to keep in mind that the latter results apply to the electron 
sector, while those established here on basis of the experimental statistical capture rate on ^He 
apply to the muon sector for couplings between the first quark generation and the second lepton 
one. On the other hand, the above lower bounds are not as stringent as those which follow 
from atomic parity violation [^], but again the latter limits apply to the electron sector and 
for neutral current interactions. It is thus fair to say that precision studies in nuclear muon 
capture on simple nuclei have at present the potential to test the SM in sectors and in ways 
complementary to those accessible through collider experiments, hence the physics interest of 
the limits established above. 

5.3 Leptoquarks 

Finally, let us consider the possibility of leptoquarks (LQ)[^, 49 1, in the notations introduced 



in Appendix 3. Among all the possible limits which may be obtained from (73), only the most 
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stringent ones are presented here. 

In the case of scalar LQ, the limit fohows from the value for the tensor coupling coefficient 
leading to 

^ > 894 GeV (95% C.L.). (94) 
A 

S VT 

The notation used here is symbolic. Indeed, referring back to the expressions for the h^^J 
coefficients in terms of the LQ couplings and parameters, one sees that the h'^_ coupling involves 
either the So(l/3) LQ with the couplings Xg*X§^, or the 5i/2(— 2/3) LQ with the couplings 
Aq* X§ . Thus in the above lower bound, M stands for the mass of one or the other of these 

'-'1/2 '-'1/2 ' 

scalar LQ, while A stands for the square root of the product of the two associated couplings 
constants. 

Similarly considering the limit in ( |73D on the coefficient h^_^_, one finds the lower bound, 

^ > 915 GeV (95% C.L.), (95) 
A 

in a similar notation referring now to either vector LQ Vo(— 2/3) or Vi/2(~l/3) (for the associated 
combination of LQ couplings corresponding to the factor A, see Appendix 3). Again, these limits 
are certainly as stringent as those recently presented in Ref . [|^, ^] . This is particularly relevant 
when one recalls again that the constraints from Ref. [p^j refer to the electron sector, while the 
limits established here apply to LQ coupling the first quark generation to the second lepton 
generation. 



6 Conclusions 

This paper provides for the first time explicit and complete analytic expressions for all obser- 
vables relevant to nuclear muon capture on a spin 1/2 isospin doublet, thus of direct use to the 
cases of the proton and ^He. The results include all possible contributions for the nuclear matrix 
elements, as well as all possible effective interactions beyond the usual electroweak charged 
interaction. The analysis was developed with great care, keeping approximations to the strictest 
minimum and in ways not affecting the final numerical results, including the calculation of 
nuclear finite size corrections to the muon overlap correction factor, and considering the limit of 
a Dirac neutrino of zero mass. Such expressions are ideally suited for tests of the Standard Model 
from precision muon capture experiments, both in addressing still open questions related to the 
strong nonperturbative sector of the quark interactions probed through the electroweak sector, 
as well as in probing for any new physics which is lying in store just waiting to be discovered 
experimentally. 

These results also enable specific predictions for observables with a precision on a par 
with the challenge on the experimental side as well. Specific new results for muon capture on 
the proton have been given, and some others in the case of ^He. This situation is particularly 
relevant given the recent precision measurement of the statistical capture rate on ^HelTl], and 



the projected 0.5% precise measurement of the singlet rate on the proton|18, 14|. 

The discussion showed how stringent limits can be inferred from the ^He result already, 
both within the SM and beyond it, sometimes in real competition with results from collider 
experiments, and often complementary to these. The physics reach of the foreseen experiment on 
the proton has also been assessed. In due time, when that experiment will have been completed 
to the desired level of precision, similar but improved limits and tests on the SM will be inferred, 
and in particular the issue of the induced pseudoscalar nucleon form factor settled once and for 
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all, presumably in perfect agreement with the theoretical expectation. Also, the potential physics 
reach of other observables may now be assessed completely on the basis of the expressions of 
this paper, even though the actual measurement of any of these observables, which all involve 
either initial polarization muonic atom states or final polarization measurements, or both, poses 
an almost impossible experimental challenge. But such a situation has never deterred any 
experimentalist at heart, quite to the contrary! We hope that this paper will also be of use in 
the experimental pursuit of the impossible polarization observables. 
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Appendix 1: The Overlap Reduction Factor 



In the calculation as outlined in Section 2, it is implicitly assumed that one is dealing with point- 
like particles whose quantum states are described by plane waves. Since this is clearly not the 
case, further corrections have to be introduced in order to account for the finite spatial extent 
of the nuclear charge distributions, for the bound state character of the initial muonic atom, 
and for possible relativistic corrections since it is the non relativistic Coulomb wave function 
ipcif) which is initially considered to represent the muon probability amplitude rather than the 
full-fledged solution to the relevant Dirac equation. Indeed, all these effects ought to be carefully 
assessed in order to obtain a trustworthy evaluation of the reduction factor C introduced in @ , 
with a precision at the required level with respect to experimental aims. 

This Appendix presents such an evaluation, following closely the discussion developed in 
Ref. in the case of ^He. Since such an analysis, which is important for our purposes, is not 
available in the literature, we felt it useful to include it here, the more so since it seems to have 
been the intention of the author of Ref.Q to make it available. 

Nuclear charge distributions, whether of electric or electroweak matter, need to be modeled 
in such an analysis. This we shall do in terms of a spherically symmetric density Pmir), where 
the index distinguishes the different types of matter encountered in the problem, normalized to 
unity over the volume of space, 

d^rpm{r) = 1, (96) 



/(<X)) 

and dependent on a single length scale parameter directly related to the mean square radius 
of the distribution. The simplest such model which we shall use is of the form 

such that 

=< r2 >„= / dVrVm(r) = 12a^. (98) 

J(oo) 

Other charge distribution models may be considered of course leading to no significant dif- 
ference in the evaluation of the reduction factor C. Note that the form factor in momentum 
space associated to the model ( |97| ) is of the usual dipolar form, 

^ra{q^) = n (99) 

(1 - ai^q^Y 

Given such models, the reduction factor Cm associated to each of these nuclear matter 
distributions is thus given byS 



V'i(r) 



d'r{l + v'ai,Ypm{r)jo{yr) (100) 

(oo) V'l (0) 

where v stands for the neutrino energy. In this expression, the factor (1 -|- v'^a^Y stems from 
the normalization of the form factor Fmiq^) or the normalization condition for the distribution 
Pm{r), jo{i^r) = sm{h'r) / (vr) is the spherical Bessel function associated to the angle-integrated 
neutrino plane wave function e~*^'^, tpi{r) is the ground state wave function of the muonic atom, 
and finally 'il^f^\r) = ipcif) is the IS ground state Coulomb wave function of the muonic atom 
since this specific choice was made in the normalization of the capture distribution (|2|). Hence, 
only the function ^l^l{r) is left to be computed in order to evaluate the above overlap integral 
giving the reduction factor Cm- 
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The wave function ^i{r) is that of the muon ("carrying" the reduced mass fi of the muonic 
atom) in the electrostatic field of the nucleus of finite size. In principle, this function is to be 
obtained by solving the associated Dirac equation in the given electrostatic potential, but one 
would expect that the ensuing relativistic corrections should be sufficiently small to be ignored 
since of order (qZ)^, which is the squared velocity of the bound muon. This expectation is 
indeed borne out by the detailed numerical resolution of the associated Dirac equation in the 
case of ^He[^ and of the proton, and this relativistic correction shall thus not be included here 
in the evaluation of Cm- 

Hence, this leaves only to solve the Schrodinger equation in the same electrostatic potential 
of the electric charge distribution of the initial nucleus. For the latter distribution, we shall 
again use a model of the form (|9^), with a parameter Oc related to the electric charge mean 
square radius of that nucleus. Correspondingly, in addition to the usual Coulomb potential for 
a point charge of value (+Ze), one needs to include in the Schrodinger equation the following 
perturbation in the potential^, 

AV{r) = aZfic (- + e~'^/"= . (101) 



r 2a, 



The contributions of this term to the IS ground state wave function ipi{r) are then computable 
through perturbation theory. It appears that the first order correction is already sufficient for 
our purposes since it is proportional to the factor (cc/ao)^ ~ 10~^, where qq is the usual atomic 
Bohr radius for the muonic atom, thus given by 



1 he , ^ 



The spherically symmetric solutions to the Coulomb problem are well known. For the nS 
state, n = 1, 2, • • •, one has 

^f(r) = ^i=e-'-/("-)4-,(^) , Vc(0)^4°^(0) = ^, (103) 
where L\_i{x) are the usual Laguerre polynomials. The associated (binding) energy eigenvalues 



are 



4°^ = -i(aZ)V'A. (104) 

Adding the correction (|101| ), first-order perturbation theory then leads to the following 15* 
ground state normalized wave function. 



V,i(r)=Vr(r)-16U Y.^,-, rV'i^^W, (105) 

Vao/ „=2Vn31-^ 

while the associated energy is 




^i = -i(aZ)Vc'fl-16(^^) ). (106) 



The evaluation of the overlap reduction factor Cm in ( |10[)|) is now straightforward using 
the integrals 



cPrpm{r)i.^\r) = ^ "° (l - ^) . (107) 

\ ' n ao J 



n~2 
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One then finds 



( 1 4- '^"i 



1 + ly'^al 



1 + 



v 



(108) 



where in the second term which follows from the first-order correction in ( |105| ) the hmit am 
has been taken. Numerically, this result does not differ from the same expression expanded to 
second order^ in Oc/ao and Om/aoQ, 



Cn 



1 



1 



3^ ""m 



1 + 1/2qi2 



^1+6 
«0 



1 



(109) 



Let us now turn to numerical evaluations, first in the case of ^He. The values for the rele- 
vant parameters may only be inferred from a nuclear model calculation of the corresponding 
mean square charge radii. This is done in Ref.P] with the following results. 



1(1) 



0.554 fm 



Ha) 



0.512 fm, 



(110) 



where a(i) (resp. a(o-)) stands for the parameter Um associated to form factors for the vector 
(resp. axial) current (thus including the electromagnetic current in the vector case). Using then 



the mass values and neutrino energy given in Section 3.1, one finds 



^ = 4.17-10-3 



(1) 



0.9777 



0.9790. 



fllll 



Weighing these two reduction factors with the relative vector and axial current contributions 
to the capture rates, an effective value of Ces = 0.9788 is obtained, thus finally leading to the 
overlap reduction factor for ^He as determined in Ref.[0], 



C(3He) = 0.979. 



fll2) 



In the case of the proton, things may be done with great precision since a great deal 



Given the value 



is known about the nucleon electromagnetic form factors] 42, ^ 
0.847 fm||4^ (known to 1% precision) for the proton electric charge distribution, one has 

^ = 8.59 • 10-^ 

Co 



(113) 



On the other hand, the values for Cv,m,a niay be determined using the associated charge radii 



and rA discussed in Section 4.1, including their uncertainties. One then finds, 



Cv = 0.99534 ± 0.00004 , Cm = 0.99460 ± 0.00005 , Ca = 0.9960 ± 0.0002. 



(114) 



Assuming then that the correction factor associated to the nuclear matter distribution characte- 
rized by the induced pseudoscalar form factor gp is also given by Ca, and weighing each of these 
correction factors with the respective relative contributions of the associated form factors to the 
capture rates, the following effective value for the overlap reduction factor is finally obtained. 



C{p) = 0.9956. 



(115) 



Note that this value agrees essentially with the 0.4% correction quoted in Ref.|21| in the case of 
hydrogen. Here however, a careful assessment of all possible effects has led to this final value, 
with increased precision. 

^°Both these ratios being on the order of 10~^, cubic corrections are indeed totaUy neghgeable for our purposes. 
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Appendix 2: Polarization and Hyperfine States 

In the calculation outlined in Sect.|2|, the natural representation of polarization states is through 
the normalized spin vectors s^^i^2 rather than the hyperfine states (5 = l,m = 0, ±1) and 
(5 = 0, m = 0) of the muonic atom. However, there exists a correspondence between these two 
bases for the initial spin degrees of freedom, and the purpose of this Appendix is to indicate how 
results pertaining to these hyperfine states may be extracted from the results derived in Sectj^. 

First, let us note that for a spin 1/2 system, any of its normalized quantum states may 
be parametrized as 

\s >= e-'"^/\os^\+ > +e'^/^sm^\->, (116) 

where |ib > represent the basis of states with spin eigenvalues m = ±1/2 with respect to 
an arbitrary quantization axis. The angular variables {6, ip) may be interpreted as being the 
spherical coordinates for a unit vector s whose components are defined as in (15), since one finds 
for the associated spin component operators given in terms of the usual Pauli matrices a, 

< a >=< s\a\s >= s (117) 

(note also that the state \s > corresponds to the bi-spinor x+{s) in (|16|)). 

Let us now consider the quantum spin states of a bound system of two spin 1/2 particles, 
such as the muonic atom, obtained through the tensor product >^ |si >i of the associated 
spin states, 

|^M>M=e-*'''^^'cos^|+>^ +e*'^-/2gin^|->^, (118) 

|si >i=e-^^i/2posy |+>i +e^^^l'^siu^-^\->i . (119) 
In terms of the usual \S,ra> hyperfine state 

|1,1 >= 1+ >^ 1+ >i , |1,-1 >= I- >^ I- >i, (120) 

|1,0 >= [|+ >^ I- >i +1- >^, 1+ >i] , |0,0 >= [|+ >^ I- >i -I- 1+ >i] , (121) 
one then finds 

\si >i=^i,i|l,l > +Ai,_i|l,-1 > +^i,o|l,0 > +Ao,o|0,0 >, (122) 

where 

^^^^ = e-^(^''+^i)/2c^ci , = e^(^^+^^)/2s^si, (123) 

^1,0 = ^ \e-''^^^-^^~^'^c^si + e^('^^-^i)/2^^^il , (124) 
V2 L J 



A. 



0,0 



1 

V2 



(125) 



with = cos^^/2, = sin 0^/2, ci = cos0i/2 and si = sin^i/2. In particular, the hyperfine 
populations are then given as 

iVi,i = J[l + cos0^][l + cos0i] , A^i,_i = ^[l-cos0^][l-cos0i], (126) 
-^1,0 ~ [■*■ ~ ^1 + cos(99^ — (/^i) sin 6^ sin ^i] , (127) 
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0,0 



[1 — cos 6^ cos 9i — cos(</?^ — ^i) sin 6f^ sin ^i] 



(128) 



with Ni^i + TVi _i + Nifi + iVo,o = 1, as it should. 

Let us now consider the matrix element M.\{Sfj,, si) =< §2, X\Hes\s^, si > which was 
expressed in Sect. |2.2| (in a notation which should be self-explanatory). In terms of the above 



change of basis in spin space, one then finds the following decomposition into hyperfine capture 
amplitudes A^^'™", 



(129) 



The hyperfine capture distributions are then given as in (Q) with the quantity |A^aP replaced by 
lA^^'^p. However, the former quantity evaluated in terms of the latter also involves interference 
contributions from different hyperfine amplitudes, which must be disposed of. 

In the case of the (S = l,m = ±1) hyperfine states, this is straightforward, since one has 



\M 



1,±1|2 
A I 



\Mxi±e3,±e3)\\ 



(130) 



where the right-handed orthonormalized basis {61,62,63} with respect to which the spherical 
coordinates {9, ip) are defined, has been introduced. 

The situation for the two other hyperfine states {S = l,m = 0) and {S = 0,m = 0) is 
more involved however, because of their intertwined character in terms of the >^ l^i >i spin 
states. One way to proceed is as follows. 

Consider the following specific combinations, 



X = \Mxie3,-e3)\'^ + |A^A(-e3,e3) 



1,0,2 ^|_^0,0|2^ 



Y 



|A^A(ei,ei)|2 - \Mxiei,-ei) 



-ei,6i 



l-^A 



A 



and 



Z 



-|A4a(cos (pfj^ei + sin 93^62, cos (^161 + sin 93162) 



\M 



1,0|2 



\M 



0,0|2 



2Re (^m]:^m]:'^ 



-ei, -61J 



\iPfj.—0,ipi=Q 



(131) 
(132) 

(133) 



Then it follows obviously that the remaining two hyperfine state capture distributions are ob- 
tained from 



i^^aT 



x + l(y + 4Z) 



i-Ma'T 



X--{Y + 4Z) 



(134) 



Applying this procedure to the quantity A/a defined in (|42|) then leads to the hyperfine 
state distributions detailed in Sect. 2^. Similarly, the quantity Vx defined in (42) would lead 
to the final nucleus polarization state associated to capture from each of the muonic hyperfine 
states, as given in (E 



Appendix 3: Left-Right Symmetric Models, Leptoquarks and 
Contact Interactions 

The purpose of this Appendix is to provide explicit expressions for the effective coupling coeffi- 
cients hj^^ introduced in (^) in terms of the parameters of specific models or parametrizations 
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for physics beyond the SM. This allows for model independent bounds that could been deter- 

S VT 

mined for the h^^' coefficients from some given experiment, to be translated into may be more 
physically tangible numbers to be compared with the reach of other experiments, especially at 
high energy colliders. Three general classes of such models beyond the SM are considered here, 
namely left-right symmetric models (LRSM), contact interactions and lepto quarks E^, 49 1. 



Left-Right symmetric models 

In the case of LRSM, we refer to the notations, discussion and references in Ref.|^^. For 
the process of interest in this paper, in the limit that Higgs exchange is ignored, which ought 
to be justified given the necessarily large masses for such particles in LRSM as well as their 
small Yukawa couplings directly proportional to the masses of the quarks and leptons involved 
in the process, only V and A couplings arise in addition to the purely {V — A) character of 
the SM charged electroweak interactions. Two massive charged gauges bosons appear, 
with masses Mj^, which, up to their mixing through an angle (" possibly accompanied by a 
CP violating phase uj, are associated to {V — A) and {V + A) interactions characterized by 
gauge coupling constants gL,R- In tbe fermionic sector, flavour mixing is also parametrized by 
Cabibbo-Kobayashi-Maskawa mixing matrices, associated to each chirality sector, denoted V^'^ 
and U^'^ for the quark and lepton sectors, respectively (indeed in a generic LRSM, neutrinos 
are massive and thus flavour- mix with one another as the quarks do). Let us introduce the 
following combinations of parameters. 

Given the parametrization @, one then finds for the only non vanishing effective coupling 
coefficients, 

2 2 

^Vud hl^ = p (l + St") Uf:, , ^Vud h^_^ = -prt (1 - 6) e'^ (136) 

= pr^ (t" + 6) C/«z;„rf , = -prt (1 - 6) e— Uf},. (137) 

Here, the indices i = 1,2,3 and p on the neutrino CKM matrix elements U^^ stand for the 
neutrino mass and muon flavour eigenstates, respectively. 

Contact interactions 

Let us now turn to the parametrization of so-called contact interactions, which provide a 
simple-minded model to probe the scale for compositeness of quarks and leptons. Such interac- 
tions are typically represented through an effective four-fermi coupling of the form|^^, say in 
the case of vector operators, 

4^,1,2 I^V^TM^^i^V^T^^r,^^', (138) 

where £.^^^,(2 = =^1) 9c is a contact interaction coupling constant, A is the associated energy scale, 
and Prii2 are the chirality projectors already introduced in (^). In the case of nuclear muon 
capture, the different spinor fields appearing in this expression correspond of course to those of 
(^), namely the u and d quarks, and the muon and its associated neutrino. It is conventional [p3|| 
to fix the scale A by setting = Att (which in the case of QED would amount to having the 
fine structure constant set to unity, a = 1). 

Clearly, such contact interactions may be introduced for any of the scalar, vector and 
tensor couplings included in (H), and for any combination of the fermion chiralities involved. 
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Thus, one could associate a scale A say to vector interactions of LL, VV, VA, etc ... chiralities, 
in any combination possible, and similarly for scalar and tensor interactions. Note that in the 
case of nuclear muon capture, these contact interactions are all related to processes which couple 
the quarks of the first generation to the leptons of the second. 



Leptoquarks 

Finally, let us consider the case of leptoquarks (LQ)[^, 33]. These particles come in 
two varieties, namely either spin or spin 1, and their name derives from the fact that they 
couple always to a quark and a lepton in a single vertex. When one allows for right-handed 
neutrinos as well (which was not considered in the original discussion |^^), there are six different 
types of scalar and of vector LQ, characterized by their weak isospin and electric charge. Their 
quantum numbers under SU{3)c x SU{2)i x U{1) are as follows, for scalar LQ, 



(139) 



and for vector LQ 
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(140) 



where each time the electric charge content of the associated isospin multiplet is given on the 
right. Note that the lower index carried by each of these fields labels its weak isospin value. 
Moreover, Sq^, and Vqi, are those LQ related to the introduction of right-handed neutrinos. 

The scalar LQ interactions are then parametrized according to the Lagrangian density, 



J^S = X^So^T2iLSl _ _ 

+ >^§o^RmSl + x§y'kmsl + Xld%^^RSl + Af: 



+ >^s^/^urS\i^Il + X§^^^qLSl^2^T2HR + X^^^^dRS\i^lL + \f^^qLSl/2ir2VR 
+ ^sMSi ■ infii + h.c. 
and similarly for vector LQ, 



(141) 



VOv 



dR-ff,I^RVi^J 



+ Xv^WVi^^ ■ If^riL + h.c. 



(142) 



Here, L and R stand for definite chiral components of the spinor fields, the upper index 
refers to the charge conjugate fields, and qi and II stand for the following quark and lepton 
doublets 



qL ■■ 



(143) 
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Finally, the different \gy coefficients are complex constant parameters, the LQ coupling con- 
stants. In general, these coefficients are matrices in generation space, but in the case of nuclear 
muon capture, only those LQ couplings between the first quark generation and the second lepton 
generation, including the possibility of right-handed neutrinos, are relevant, hence our choice of 
notation. 

S V T 

Given these definitions, the induced non vanishing effective couplings coefficients in 
(^) are expressed as, 

SM'^'^'y") 2M2J-2/3) 2M2^^^(-1/3)' ^ ' 

2 T \L \Ru* 1 X^* 

8M^ y -+) ~ + 8 (1/3) ^ 8M|^^(-2/3)' ^ ^ 

8M2 '^"'^ r+-J - 8M|^(l/3) 8M|^^(-2/3)' ^ ' 



-.2 



1 ^uAv!^* 1 ^Vi/->^Vt/- 



Vud[n^+) - 2M2f-2/3l 2 Ml f-l/3V ^ 



8^2 V ++y 2 Ml (-2/3) 2M? (-1/3) 

" ''1/2 



2 2 1 I \L |2 1 I \L |2 1 I \L |2 i | \L |2 

9 y (.V Y ^ 9l ySM , 1 l^gpl _ 1 1^1 , 1 I Vol _ 1 IVll 

8M2 ""^ V — y 8M2, ""^ ^8M2^(l/3) 8M|^(-l/3) 4M2^(-2/3) 4M2^(-2/3)' 

(148) 

„2 1 \^ \Rv* -I \R \Ru* 

8M2 '^"'^ 8 M2^(1/3) ^ 4 M4(-2/3) ' ^ ' 

2 1 1 \i' x-Ri^* 1 A^* 

^ K.ifc)- = -i,4?^ + ^.#^. (150) 



8M2 "'='2V-+/ 32M|^(1/3) 32M|^^(-2/3)' 
1 AT Y_ , 1 A^^Ag, 1 X's:,,Xl, 



In these expressions, which of the LQ isospin component contributes to a each coefficient is 
indicated by giving its electric charge in the parenthesis following the mass value in the mass 
contributions l/Mgy. Note that the Sqi, and Vq^ LQ do not contribute to these coefficients. 
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a(AS;Fx) 






ct(Aa;Fx) 


<y{A^;Fx) 


(T{At-Fx) 


Fv 


-0.726 


+0.790 


+0.301 


-4.572 


+0.798 


+0.0616 


Fm 


+0.420 


+0.233 


+0.293 


+0.561 


-0.287 


+0.254 


Fa 


+2.621 


+0.997 


+1.521 


+4.894 


-0.134 


-1.063 


Fp 


-0.314 


-0.0210 


-0.116 


-0.884 


-0.377 


+0.747 


Fs 


-0.0155 


+0.0178 


+0.0071 


-0.101 


+0.0173 


+0.0017 


Prp 


-0.0157 


-0.00112 


-0.0058 


-0.044 


-0.0185 


+0.0368 



Table 1: Sensitivities of the different observables O = y^^''^'^^^^ , A^, A^, At for muon capture on 
^He, with respect to the nuclear form factors associated to the vector and axial quark operators. 
See (66) for the definition of a{0;Fx)- 





ct(AS;/ix) 


a{\^;hx) 


a(A^*^*;/ix) 


a{A/^]hx) 


(j{Ay; hx) 


a{Au hx) 




-0.840 


+0.946 


+0.370 


-5.38 


+0.901 


+0.126 




-0.810 


+0.948 


+0.381 


-5.30 


+0.937 


+0.0544 


hi 


-0.825 


+0.947 


+0.375 


-5.34 


+0.919 


+0.0902 


hi 


-0.0151 


-0.00093 


-0.00563 


-0.0426 


-0.0180 


+0.0361 




+10.1 


+3.78 


+5.82 


+19.00 


-0.733 


-3.83 


h^ 


+2.00 


+2.00 


+2.00 


+0.00 


+0.00 


+0.00 


/.^ 


-2.61 


+0.0478 


-0.810 


-8.022 


+1.022 


+0.632 



Table 2: Sensitivities of the different observables O = y?'"^'^^^^ , A^, A^^ At for muon capture on 
^He, with respect to the effective coupling coefficients ^. Sensitivities for those couplings 
which do not appear in this table are identically vanishing. The scalar and pseudoscalar com- 
binations are defined by hi = hl_ + /i^^ and hi = hl_ — while it is understood that all 
form factors Gs, Gp and Gt are set to unity to obtain the numbers in this table. See (36) for 
the definition of £7(0; hx)- 
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o-(^a;5x) 


a{A^;gx) 


(^{At;9x) 


9V 


+0.466 


-1.129 


+0.386 


+0.108 


+13.16 


-2.39 


9M 


+0.151 


+0.680 


+0.177 


-0.0357 


-0.317 


+0.551 


9A 


+1.567 


+1.440 


+1.561 


+0.00856 


-18.85 


+0.991 


9P 


-0.184 


+1.008 


-0.125 


-0.0804 


+6.01 


+0.844 


9S 


+0.0232 


-0.0718 


+0.0185 


-0.00641 


+0.724 


-0.139 


9t 


+0.0238 


-0.125 


+0.0164 


+0.01003 


-0.759 


-0.105 



Table 3: Sensitivities of the different observables O = X^''^'^^'^^ , A/^, Ay, At for muon capture on 
the proton, with respect to the nuclear form factors associated to the vector and axial quark 
operators. See (|6^) for the definition of a{0;Fx)- 





ct(AS;/ix) 


cT{X^;hx) 




a{A/^;hx) 


a{Ay; hx) 


a{Au hx) 




+0.434 


-1.394 


+0.343 


+0.123 


+12.16 


-2.562 




+0.391 


-1.157 


+0.314 


+0.104 


+13.57 


-2.363 


hi 


+0.412 


-1.274 


+0.328 


+0.114 


+12.86 


-2.461 


hi 


+0.0217 


-0.118 


+0.0147 


+0.00945 


-0.710 


-0.0986 




-5.21 


-5.58 


-5.23 


+0.0251 


+55.21 


-3.65 




+2.00 


+2.00 


+2.00 


+0.00 


+0.00 


+0.00 




-0.767 


-2.90 


-0.872 


+0.144 


+25.68 


-3.669 



Table 4: Sensitivities of the different observables O = }?''^'^'^^^ , A^, A^, At for muon capture 

S V T 

on the proton, with respect to the effective coupling coefficients hj^^ . Sensitivities for those 
couplings which do not appear in this table are identically vanishing. The scalar and pseudoscalar 
combinations are defined by hi = hl_ + hl_^_ and hi = hl_ — while it is understood that 
all form factors Gs, Gp and Gt are set to unity to obtain the numbers in this table. See ( |66[ ) 
for the definition of a{0; hx)- 
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